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ABSTRACT 



Context. In a strongly stratified turbulent layer, a uniform horizontal magnetic field can become unstable to form spontaneously local 
flux concentrations due to a negative contribution of turbulence to the large-scale (mean-field) magnetic pressure. This mechanism is 
of interest in connection with dynamo scenarios in which most of the magnetic field resides in the bulk of the convection zone, and 
not at the bottom, as is usually assumed. Recent work using the mean-field hydromagnetic equations has shown that this negative 
effective magnetic pressure instability (NEMPI) becomes suppressed at rather low rotation rates with Coriolis numbers as low as 0.1. 
Aims. Here we extend these earlier investigations by studying the effects of rotation both on the development of NEMPI and on the 
effective magnetic pressure (turbulent and non-turbulent contributions). We quantify the kinetic helicity resulting from rotation and 
stratification and compare with earlier work at smaller scale-separation ratios. We also determine the sensitivity of surface diagnostics 
of magnetic helicity. 

Methods. We use direct numerical simulations (DNS) and mean-field calculations of the three-dimensional hydromagnetic equations 
in a Cartesian domain and analytical studies using the r approach. 

Results. We find that the growth rates of NEMPI in earlier mean-field calculations are well reproduced with DNS and that the 
rotational effect on the effective magnetic pressure is negligible as long as the production of flux concentrations is not inhibited by 
rotation. In that case, kinetic and magnetic helicity are also found to be weak. 

Conclusions. Production of magnetic flux concentrations through the suppression of turbulent pressure appears to be possible only 
in the upper-most layers of the Sun, where the convective turnover time is less than 2 hours. 

Key words, magnetohydrodynamics (MHD) - hydrodynamics - turbulence - Sun: dynamo 



1. Introduction 

In the Sun, magnetic fields are produced by a large-scale dy- 
namo (see, e.g., Moffatt, 1978; Parker, 1979; Krause & Radler, 
1980; Zeldovich et al., 1983; Ossendrijver, 2003; Brandenburg 
& Subramanian, 2005a). Although many details of this process 
remain subject to debate, it seems relatively clear that rotation 
enhances the efficiency of the dynamo if the Coriolis parame- 
ter is not very large. In the absence of rotation and shear, only 
small-scale magnetic fields are generated by what is often re- 
ferred to as small-scale dynamo action (see, e.g., Zeldovich et 
al., 1990; Brandenburg & Subramanian, 2005a). Rotation leads 
to an a effect (Steenbeck et al., 1966) if there is also stratifica- 
tion in density or turbulent intensity. The a effect can produce 
mean magnetic field and net magnetic flux. 

Stratification leads to yet another effect which does not 
produce magnetic flux, but merely concentrates it locally by 
what is now referred to as negative effective magnetic pressure 
instability (NEMPI). Direct numerical simulations (DNS) of 
Brandenburg et al. (201 la) have shown in surprising detail many 
aspects of NEMPI that were previously seen in mean-field sim- 
ulations (MFS) of Brandenburg et al. (2010) and that have been 
anticipated based on analytical studies for some time (Kleeorin 
et al., 1989, 1990, 1993, 1996; Kleeorin & Rogachevskii, 1994; 
Rogachevskii & Kleeorin, 2007). 

The main physics of this effect is connected with the sup- 
pression of turbulent pressure by a weak mean magnetic field 



that is less than the equipartition field. At large Reynolds num- 
bers, the resulting reduction of the turbulent pressure is larger 
than the added magnetic pressure from the mean magnetic field 
itself, so that the effective magnetic pressure that accounts for 
turbulent and non-turbulent contributions, becomes negative. In 
a strongly stratified layer, i.e., a layer in which the density varies 
much more rapidly with height than the magnetic field, this 
leads to an instability that is analogous to Parker's magnetic 
buoyancy instability, except that there the magnetic field varies 
more rapidly with height than the density. Because the effective 
magnetic pressure is negative, magnetic structures are negatively 
buoyant and sink, which has been seen in DNS of Brandenburg 
etal. (2011a). 

One of the main successes of recent comparative work be- 
tween DNS and MFS is the demonstration of a high degree of 
predictive power of MFS. The examples include details regard- 
ing the shape and evolution of structures, the dependence of their 
depth on the magnetic field strength, and the dependence of the 
growth rate on the scale separation ratio. Recent MFS of Losada 
et al. (2012) (hereafter LBKMR) have shown that in the pres- 
ence of even just weak rotation, the growth rate of NEMPI is 
significantly reduced. Expressed in terms of the Coriolis num- 
ber, Co = 2f2/u rms fcf, where il is the angular velocity, u lms is 
the rms velocity of the turbulence, and fcf is the wavenumber of 
the energy-carrying eddies, the critical value of Co was predicted 
to be as low as 0.03. Although this value does not preclude the 
operation of NEMPI in the upper parts of the Sun, where Co 
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is indeed small (about 10~ 4 at the surface), it does seem sur- 
prisingly low, which raises questions regarding the accuracy of 
MFS in this case. The purpose of the present paper is therefore 
to compare MFS of LBKMR with DNS of the same setup. 

The MFS of LBKMR were done for a case for which com- 
parisons between MFS and DNS do exist in the special case of 
Co = (Kernel et al., 2012a). The aim of this paper is to extend 
these effects to the case of non-vanishing values of the Coriolis 
number Co. We use DNS of the MHD equations and analytical 
studies based on the r approach. We recall that we adopt here 
an isothermal stratification and an isothermal equation of state. 
This is done because the effect that we are interested in exists 
even in this simplest case where temperature and pressure scale 
height are constant. Non-isothermal setups have been studied at 
the mean-field level both with (Kapyla et al., 2012, 2013) and 
without (Brandenburg et al., 2010) entropy evolution included. 
In a stably stratified layer, entropy evolution leads to an addi- 
tional restoring force and hence to internal gravity waves (Brunt- 
Vaisalit oscillations) that stabilize NEMPI (Kapyla et al., 2012). 
Thus, by using both isothermal stratification and an isothermal 
equation of state, we recover a situation that is similar to an adi- 
abatic layer, except that then the temperature and hence the pres- 
sure scale height decrease with height. 

2. The model 

We consider DNS of an isothermally stratified layer 
(Brandenburg et al., 2011a; Kernel et al., 2012a) and solve 
the equations for the velocity U, the magnetic vector potential 
A, and the density p, in the presence of rotation il, 



— = -20xt/-^Vlnp+-JxB + /+g + F„,(l) 



Df7 

dA 

~dt 
dp 

dt 



U xB + r]V 2 A, 
-V • pU, 



(2) 
(3) 



where D/Di = d/dt + U ■ V is the advective derivative, v 
is the kinematic viscosity, rj is the magnetic diffusivity due to 
Spitzer conductivity of the plasma, B = Bq + V x A is 
the magnetic field, B = (0, B ,Q) is the imposed uniform 
field, J = V x B/po is the current density, p is the vac- 
uum permeability, F v = V • (2vpS) is the viscous force, Sjj — 
^(djUi + diUj) — ^Sij'V -U is the traceless rate-of-strain tensor. 
The angular velocity vector ft is quantified by its scalar ampli- 
tude ft and colatitude 9, such that fi = £1 (— sin 9, 0, cos 9). As 
in LBKMR, z corresponds to radius, x to colatitude, and y to az- 
imuth. The forcing function / consists of random, white-in-time, 
plane, non-polarized waves with a certain average wavenumber 
kf. The turbulent rms velocity is approximately independent of 
z with u rms = (it 2 ) 1 / 2 k, 0.1 c s . The gravitational acceleration 
g = (0,0, —g) is chosen such that k\H p = 1, so the density 
contrast between bottom and top is exp(27r) s=s 535 in a domain 
— 7r < k\z < 7r. Here, H p = c 2 /g is the density scale height 
and fci = 2tt/L is the smallest wavenumber that fits into the cu- 
bic domain of size L 3 . In most of our calculations, structures 
develop whose horizontal wavenumber k x is close to k\. We 
adopt Cartesian coordinates (x,y,z), with periodic boundary 
conditions in the x- and y-directions and stress-free, perfectly 
conducting boundaries at the top and bottom (z = ±L z /2). In 
all cases, we use a scale separation ratio kf/ki of 30, a fluid 
Reynolds number Re = u lms / vk{ of 36, and a magnetic Prandtl 
number PrM = vjr\ of 0.5. The magnetic Reynolds number is 



therefore ReM = Pr^/Re = 18. The value of B is specified in 
units of the volume-averaged value B cq0 = y/poPo u ms , where 
po = (p) is the volume-averaged density, which is constant in 
time. In our units, k\ = c s = po = po = 1. In addition 
to visualizations of the actual magnetic field, we also monitor 
B v , which is an average over y and a certain time interval At. 
Time is sometimes specified in terms of turbulent-diffusive times 
tr]tokj, where rj t o — u rms /3kf is the estimated turbulent diffu- 
sivity. 

The simulations are performed with the Pencil Code 
http://pencil-code.googlecode.com which uses 
sixth-order explicit finite differences in space and a third-order 
accurate time stepping method. We use a numerical resolution 
of 256 3 mesh points. 

We compare with and extend earlier MFS of LBKMR, where 
we solve the evolution equations for mean velocity U, mean 
density ~p~, and mean vector potential A, in the form 



dU 

~dt_ 

dA 

~dt 

&p 

dt 



= —U ■ VU - 2CI x U - c 2 s V lnp + g + JF MK , (4) 



UxB-{r]t + Tf)J, 
-U-Vp-pV ■ U, 



where J M k = + ^k. with 

being the mean-field magnetic pressure force, and 
= {vt + v) (V 2 U + i VV -U + 2SV lap) 



(5) 
(6) 

(7) 
(8) 



is the total (turbulent plus microscopic) viscous force. Here, 
— + — ^SijV ■ U is the traceless rate-of-strain 

tensor of the mean flow and q p is approximated by (Kernel et al., 
2012b) 



PI 



(9) 



which is only a function of the ratio (3 = \B\/B cq . Here, /3* 
and /3 p are coefficients that have been determined from previous 
numerical simulations in the absence of rotation (Brandenburg et 
al., 2012a). In Eq. (7) we have taken into account that the mean 
magnetic field is independent of y, so the mean magnetic tension 
vanishes. 

The strength of gravitational stratification is characterized 
by the nondimensional parameter Gr = g/c 2 k{ = (Hpkf)^ 1 
(Brandenburg et al., 2012b). Another important nondimen- 
sional parameter is the Coriolis number, Co = 2f2/u rms fcf. 
Alternatively, we normalize the growth rate of the instability by 
a quantity 



(10) 



which is motivated by the analytic results of LBKMR and the 
finding that NEMPI is suppressed when 217 ^ A»o- 

3. Coriolis effects of NEMPI in DNS 

We have performed DNS for different values of Co and calcu- 
lated the growth rate A; see Fig. 1, where we compare with the 
MFS of LBKMR, who used q p0 = 20 and /3 p = 0.167. This 
set of parameters is based on a fit by (Kernel et al., 2012a) for 
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magnetic field. This allows us to compute q p (f3) using Eq. (17) 
of Brandenburg et al. (2012a): 



Fig. 1. Dependence of A /A,o on 2tt/ A, o for DNS, compared with MFS 
(i) where g p o = 20 and /3 P = 0.167 and MFS (ii) where q p o = 32 and 
P p — 0.058. In this case no growth was found for Co > 0.03. In all 
cases we have Bo /B cci o = 0.05. 



kf/ki = 30. Note that the growth rates for the MFS are about 3 
times larger than those of the DNS; see Appendix A for a revised 
comparison of the 2-D and 3-D data of LBKMR. This might be 
caused by an inaccurate estimate of the mean-field coefficients. 
Indeed, according to Brandenburg et al. (2012a), these parame- 
ters should be q p o = 32 and (3 p = 0.058 for Re m = 18. We have 
now performed new 3-D MFS with this set of parameters. Those 
results are also shown in Fig. 1. It turns out that with these new 
parameters the resulting growth rates are indeed much closer to 
those of the DNS, suggesting that the former set of mean-field 
coefficients might indeed have been inaccurate. A reason for this 
might be the fact that for k[ jk\ = 30 NEMPI is very strong and 
leads to inhomogeneous magnetic fields for which the usual de- 
termination of mean-field coefficients, as used by Brandenburg 
et al. (2012a), is no longer valid. 

In Fig. 2 we show the yt-averaged B y for Co = 0.006, 
0.03, and 0.06 at different times. These figures show the gen- 
eration of structures that begin to sink subsequently. However, 
for Co = 0.03 and larger, this sinking is much less prominent. 
Instead, the structures remain confined to the surface layers, 
which is seen more clearly in visualizations of B y at the pe- 
riphery of the computational domain for Co = 0.03; see Fig. 3. 

To our surprise, the large-scale structures remain still inde- 
pendent of the r/-direction, which is clearly at variance with re- 
sults of the corresponding MFS. Even at other angles such as 
8 = 45° and 90°, no variation in the y-direction is seen; see 
Fig. 4. The reason for this discrepancy with the corresponding 
MFS is not yet understood. 



4. Effective magnetic pressure 

In this section we study the effect of rotation on the function 
q P (/3). We consider first the results of DNS and turn then to an 
analytical treatment. 



4. 1. Numerical results 

In the MFS of LBKMR we assumed that V e fi(0) does not 
change significantly with Co in the range considered. With DNS 
we can compute VcS (/?) by calculating the combined Reynolds 
and Maxwell stress for a run with and one without an imposed 



q p = -2 p(u 2 x -ul x ) + \b 2 -b% 



B 



(11) 



where the subscripts indicate values obtained from a reference 
run with Bo = 0. This expression does not take into account 
small-scale dynamo action which can produce finite background 
magnetic fluctuations bo. The effective magnetic pressure is then 
determined using the equation V c g (/?) = \[1 — q p (/?)]. The re- 
sult is plotted in Fig. 5 for three values of Co during an early 
time interval when structure formation is still weak enough so as 
not to affect the result. We also show the results using Eq. (9) for 
different combinations of /3* and f3 p . It turns out that the curves 
for Co = and 0.03 are close together, and that for Co = 0.13 
the negative pressure effect is perhaps even slightly enhanced. In 
any case, it is clear that within the range of values of Co consid- 
ered, the effect of rotation on q p is quite small. 

4.2. Theoretical predictions 

Let us now compare with theoretical predictions for q p (f3). We 
take into account the feedback of the magnetic field on the tur- 
bulent fluid flow. We use a mean-field approach whereby ve- 
locity, pressure and magnetic field are separated into mean and 
fluctuating parts. We also assume vanishing mean motion. The 
strategy of our analytic derivation is to determine the ft depen- 
dencies of the second moments for the velocity Ui (t, x) Uj (t, x), 
the magnetic field bi(t, x) bj(t, x), and the cross-helicity tensor 

bi(t, x) Uj(t, x), where b are fluctuations of magnetic field pro- 
duced by tangling of the large-scale field. To this end we use the 
equations for fluctuations of velocity and magnetic field in rotat- 
ing turbulence, which are obtained by subtracting equations for 
the mean fields from the corresponding equations for the actual 
(mean plus fluctuating) fields. 



4.2.1. Governing equations 

The equations for the fluctuations of velocity and magnetic fields 
are given by 



du(x, t) 
di 



~ + =\(b- V)B + (B ■ V)b] + 2u x fl 

P P 



(B-V)u-{u-V)B+Af b , 



(12) 
(13) 



db(x,t) 
dt 

where Eq. (12) is written in a reference frame rotating with con- 
stant angular velocity J~2, p = p' + (B ■ b) are the fluctuations of 
total pressure, p' are the fluctuations of fluid pressure, B is the 
mean magnetic field, and p is the mean fluid density. For simplic- 
ity we neglect effects of compressibility. The terms N u and TV , 
which include nonlinear and molecular viscous and dissipative 
terms, are given by 

N u = (u ■ V)u - (u ■ V)u + \ (j x 6 - j x b) + fJu), 

P 



JV = V x [u x b - u x b - ?/ V x b) , 

where p/i, (it) is the molecular viscous force and j = Vx6/^i 
is the fluctuating current density. To eliminate the pressure term 
from the equation of motion (12) we calculate Vx(Vxit). 
Then we rewrite the obtained equation and Eq. (13) in Fourier 
space. 
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Fig. 2. yt-averaged B y for Co = 0.006 (left), 0.03 (middle), and 0.06 (right) at different times. 



4.2.2. Two-scale approach 



4.2.3. Equations for the second moments 



We apply the two-scale approach and express two-point correla- 
tion functions in the following form 

Ui(x)uj(y) — J dk 1 dk 2 u i (ki)uj(k2)exp{i(ki-x + k 2 -y)} 



J dk dK fij (k, K) exp(ik-r + iK R) 
J dk (k, R) exp(ife-r) 



(14) 



(see, e.g., Roberts & Soward, 1975). Here and elsewhere, we 
drop the common argument t in the correlation functions, 

fij(k, R) = L(ui\ Uj), where 



L(a; c) = J a(k + K/2)c(-k + K/2) exp (iK-R) dK, 

with the new variables R = (x + y)/2, r = x — y, K = ki + 
ki, k = (fei —ki)j2. The variables R and K correspond to the 
large scales, while r and k correspond to the small scales. This 
implies that we have assumed that there exists a separation of 
scales, i.e., the turbulent forcing scale if is much smaller than the 
characteristic scale Lb of inhomogeneity of the mean magnetic 
field. 



We derive equations for the following correlation functions: 

fij(k,R) = L(ui;uj), hij(k,R) = Lfabj) and 
gij(k, R) — L(bi, Uj). The equations for these correlation func- 
tions are given by 



dt 
dhjj(k) 

dt 
dgij(k) 



i(k-B)^ + L? jmn f mn + J| +Af f , 



-i(k-B)^ + i* 



;)/ -i(fc-B)[/y(fc)-fttf(fc)-ft y 
+D? m (k 2 )g im (k) + I? j +rt9, 
where 

$„(fc) = p- 1 [ gij (k) - gji(-k)], 

F) n (h) — 9p- ■ O h 



L 



ijnm — D im (k\) Sjn + Dj n {k2) 5i m . 



(15) 



(16) 



(17) 



Hereafter we have omitted the U-argument in the correlation 
functions and neglected terms ~ 0(V|j), and £y„ is the fully 
antisymmetric Levi-Civita tensor. In Eqs. (15)— (17), the terms 
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Fig. 3. By at the periphery of the computational domain for Co = 0.006 (left), 0.03 (middle), and 0.06 (right) at the same times as in Fig. 2. Note 
the strong surface effect for Co = 0.03 in the last time frame. 



J\[f , Af h and M 9 are determined by the third moments appear- 
ing due to the nonlinear terms, the source terms if , Ij 1 - and 
/?• which contain the large-scale spatial derivatives of the mean 
magnetic and velocity fields, are given by Eqs. (A3)-(A6) in 
Rogachevskii & Kleeorin (2004). These terms determine turbu- 
lent magnetic diffusion and effects of nonuniform mean velocity 
on the mean electromotive force. 

For the derivation of Eqs. (15)— (17) we use an approach that 
is similar to that applied in Rogachevskii & Kleeorin (2004). We 
take into account that in Eq. (17) the terms with tensors that are 
symmetric in i and j do not contribute to the mean electromotive 
force because £ m = e m ji . We split all tensors into nonheli- 

cal, hij, and helical, h\j , parts. The helical part of the tensor 

of magnetic fluctuations h\j depends on the magnetic helicity, 

and the equation for h\ - ' follows from magnetic helicity con- 



servation arguments (see, e.g., Kleeorin & Rogachevskii, 1999; 
Brandenburg & Subramanian, 2005a, and references therein). 

4.2.4. r-approach 

The second-moment equations (15)— (17) include the first-order 
spatial differential operators applied to the third-order moments 
Af( m ). To close the system, we express the set of the third-order 
terms Sf M = M~M^ im through the lower moments M^ IV > . We 
use the spectral r approximation which postulates that the de- 
viations of the third-moment terms, MM"™ (fe), from the con- 
tributions to these terms afforded by the background turbulence, 
M M" u,0 > (k), are expressed through similar deviations of the 
second moments: 

-*M<™°>(*) = -"'■■'<*> -f"'°'< fe > 

T(k) 



5 



I. R. Losada et al.: Competing effects in concentrating magnetic flux 



t/T = 1.0QC 





B y/ B ^ 
no.3 



0.0 



-0.3 



t/T v = 1-000 ; . 




B y/ B B 

no.: 



Ifll IS 




-0.3 



t/T v = 1.303- 



'J B z% 
no.3 



f/7 1 = 1.60J 



: I - 



0.0 



-0.3 




-0.3 



no.3 




0.3 



0.3 



Fig. 4. By at the periphery of the computational domain for Co = 0.006 and 9 — 45° (upper row) and 90° (lower row), at three different times 
(from left to right). 
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Fig. 5. Normalized effective magnetic pressure, Vqs{B), for three val- 
ues of Co, compared with Eq. (9) for different combinations of 8+ and 
P P , as discussed in the text. 



(18) 

(Orszag, 1970; Pouquet et al., 1976; Kleeorin et al., 1990; 
Rogachevskii & Kleeorin, 2004), where r(k) is the scale- 
dependent relaxation time, which can be identified with the cor- 
relation time of the turbulent velocity field for large Reynolds 
numbers. The quantities with the superscript (0) correspond to 
the background turbulence (see below). We apply the spectral r 
approximation only for the nonhelical part hij of the tensor of 
magnetic fluctuations. A justification for the r approximation in 



different situations has been offered through numerical simula- 
tions and analytical studies (see, e.g., Brandenburg et al., 2004; 
Brandenburg & Subramanian, 2005b, 2007; Rogachevskii et al., 
2011). 

4.2.5. Solution of equations for the second moments 

We solve Eqs. (15)— (17) neglecting the sources 7^, If- with 
the large-scale spatial derivatives. The terms with the large-scale 
spatial derivatives which determine the turbulent magnetic dif- 
fusion, can be taken into account by perturbations. We subtract 
from Eqs. (15)— (17) the corresponding equations written for the 
background turbulence, use the spectral r approximation. We as- 
sume that the characteristic time of variation of the second mo- 
ments is substantially larger than the correlation time r(fc) for 
all turbulence scales. This allows us to get a stationary solution 
for the equations for the second-order moments, M^ u \ Thus, we 
arrive to the following steady-state solution of Eqs. (15)— (17): 



hij(k) 

9ij (k) = ir(k-B)D^ [f mj (k) - h rn3 {k)} 



-iT{k-B)<$> i3 {k), 
-1 



(19) 

(20) 
(21) 



We have assumed that there is no small-scale dynamo in the 
background turbulence. Here the operator D^- is the inverse of 
the operator Sij — tD^ (Radler et al., 2003) and the operator 

Lijmn is the inverse of the operator 5 im 5j n - r L% jmn (Elperin 
et al., 2005). These operators are given by 
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k 



(22) 



-\-J3x, Sip m £jq n kpq -\~ Bq {^imp^jpn ^jnp^-zpm)] 

o& ion kr>a ) + 0(^ 3 ), (23) 



^^imp^jqn ^pq ) 



where fc, = fej/ft, xO) = !/(! + ^ )- ^ = 2r ( fc ) ( k ' 

Si = 1 + x(2V), S 2 = Si + 2 - 4 X (^), B 3 = 2^x(2^), 

S 4 = 2x W - Si, S 5 = 2 - Si and S 6 = 2^ [x(V0 - xM, 



«7 \™) — $ij 



kikj/k , 5ij is the Kronecker tensor. 



We use the following model for the homogeneous and 
isotropic background turbulence: /|?'(fe) = (u 2 ) Pij(k) W(k), 
where W(k) = S(fc)/87rfc 2 , the energy spectrum is E(k) 
a — l)k ~ 1 (k/ko)~ q , ko — l/£f and the length if is the maxi- 
mum scale of turbulent motions. The turbulent correlation time 
is r(k) — C To (fc/fco) _M , where the coefficient C = (q — 
1 + n)/(q — 1). This value of the coefficient C corresponds to 
the standard form of the turbulent diffusion coefficient in the 
isotropic case, i.e., rj T = (u 2 ) J r(k) E(k) dk — tq(u 2 }/3. 

Here the time To = lf/y/ (u 2 ) and J (u 2 ) is the characteris- 
tic turbulent velocity in the scale £{. For the Kolmogorov's type 
background turbulence (i.e., for a turbulence with a constant en- 
ergy flux over the spectrum), the exponent \i = q — 1 and the 
coefficient C = 2. In the case of a turbulence with a scale- 
independent correlation time, the exponent /i = and the coeffi- 
cient C = 1. Motions in the background turbulence are assumed 
to be non-helical. 

Equations (19)— (23) yield: 



f lJ (k) = fg\k)-h lJ (k), 



M fe ) = 



1 + 2* 



1- V 



2(1 + 2*) 



(24) 



(25) 



where * = 2T 2 (k-c A f 
account that L^ mn P mn 



', ca = B/i/p, and we have taken into 
(fe) = Pij (k). After the integration in k 
space we obtain the magnetic tensor hij in physical space: 



(26) 



where (3 — S/S cq , and the functions <?i(/3) and qi^S) are given 
in Appendix B. We consider the case in which the angular veloc- 
ity is perpendicular to the mean magnetic field. The results can 
easily be generalized to the case of the arbitrary angle between 
the angular velocity and the mean magnetic field. 

The contribution of turbulence to the mean-field magnetic 
pressure is given by the function q p ((3) = [qi{fi) — 92 0$)] I0 2 '- 



1 



12^ 2 



.(0) 



(0), 



-2(ttT ) 2 (4 2) (0) -4C^ 2) (0) - 104 2) (4/3) 



-40Cf ) (4 i 9) 



4Ci(16/3 2 )] 



(27) 



where the functions A^'(x), C 4 u; (x), Ai(y) and Ci(y), and 
their asymptotics are given in Appendix B. Following earlier 
work (Brandenburg et al., 2012a), we now define a magnetic 
Reynolds number based on the scale if — 2%/kf, which is 
related to the Rej\/ defined earlier via Rm = 2%ReM- For 
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Fig. 6. Relative kinetic helicity spectrum as a function of GrCo 
(red symbols) compared with results from earlier simulations of 
Brandenburg et al. (2012b). 



S <C S oq /4Rm 1 / 4 , the function q p (f3) is given by 



gp (/3) = -lnRm--Co 2 , 



(28) 



and for S cq /4Rm 1/4 < S < S cq /4 the function q p (/3) is 
given by 

16 8 

<h(P) = 25 (l + 5|ln(4/3)|+ 32/3 2 )--Co 2 , (29) 

where Co = 2ilro. This shows that for the values of Co of inter- 
est (Co < 0.06), the correction to q p is negligible (below 10 -3 ), 
which is in agreement with the numerical findings in Fig. 5. 



5. Kinetic and magnetic helicity 

5.1. Helicity production 

In turbulence, the presence of rotation and stratification gives 
rise to kinetic helicity and an a effect (Krause & Radler, 1980; 
Brandenburg et al., 2013). As a measure of kinetic helicity, we 
determine the normalized helicity 



6f = IjJ ■ u/kfU 



2 

rms" 



(30) 



In Fig. 6 we compare two of our present runs at fcf/fci = 30 
with those of Brandenburg et al. (2012b) at kf/ki = 5 showing 
ef versus GrCo. For our present runs (red open symbols), kinetic 
helicity is clearly very small, which is a consequence of the small 
value of Co. Compared with earlier runs at kf/ki — 5, which 
gave £{ w 2 GrCo, the present ones show about twice as much 
helicity. Interestingly, in the limit of rapid rotation the relative 
kinetic helicity declines again when the product GrCo is larger 
than about 0.5. The maximum value of ef that can be reached is 
about 0.3. 

5.2. Surface diagnostics 

As a consequence of the production of kinetic helicity, the mag- 
netic field should also be helical. However, since magnetic he- 
licity is conserved, and since it was zero initially, it should re- 
main zero - at least on a dynamical time scale (Berger, 1984). 
This condition can be obeyed if the magnetic field is bi-helical, 
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0.010 
0.005 




-0.010 



ingful conclusions can be drawn from our results. By compari- 
son, a bi-helicity is expected in the northern hemisphere of the 
Sun, where magnetic helicity is negative on all scales except the 
largest ones where the a effect operates. In this connection we 
remind the reader that this sense is reversed in the solar wind far 
from the Sun (Brandenburg et al., 2011b). This has also been 
seen in simulations of magnetic ejecta from a dynamo-active 
sphere (Warnecke et al., 2011), which may be explained by a 
diffusive magnetic helicity transport (Warnecke et al., 2012). 
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Fig. 7. Normalized magnetic helicity spectrum for Co = 0.03 (upper 
panel) and Co = 0.13 (lower panel). 



i.e., with opposite signs of magnetic helicity at large and small 
wavenumbers (Seehafer, 1996; Ji, 1999). We may now ask 
whether signatures of this could in principle be detected at the 
solar surface. To address this question, we use our simulation 
at intermediate rotation speed with Co = 0.03, where magnetic 
flux concentrations are well developed at the surface of the do- 
main, and compare with a larger value of 0.13. 

Measuring magnetic helicity is notoriously difficult, because 
it involves the magnetic vector potential which is gauge de- 
pendent. However, under the assumption of homogeneity and 
isotropy, the Fourier transform of the magnetic correlation ten- 
sor is 



Mij(k) = (Sij - kikj) 



HoE M (k) 
ink 2 



i k k H M (k) 
8nk 2 ' 



(31) 



where k = k/k is the unit vector of k, and Em (k) and Hm (k) 
are magnetic energy and magnetic helicity spectra, which obey 
the realizability condition, 2^ ^M(fc) > k\HM(k)\, where the 
factor 2 in front of Eyi{k) is just a consequence of the factor 1/2 
in the definition of energy. Matthaeus et al. (1982) used the solar 
wind data from Voyager II io determine HM(k) from the in situ 
measurements of B and Brandenburg et al. (201 lb) applied it to 
measuring H^ik) at high heliographic latitudes where Hyi{k) 
is finite and turned out to be bi-helical. We now adopt the same 
method using Fourier transforms in the y direction. In the Sun, 
this corresponds to measuring B along a 2ir ring at fixed polar 
latitude. In Fig. 7 we show the result for two values of Co. 

It turns out that Hu{k) is compatible with zero for our in- 
termediate value of Co. For faster rotation (Co = 0.13), ifjvi(fc) 
is negative both at large wavenumbers (k 3> k{), and positive 
(but still compatible with zero within error bars) at intermedi- 
ate wavenumbers (0.15 < k/kf < 0.6). For k/kf < 0.15, 
the magnetic helicity is again negative. However, the error bars 
are large and rotation is already so fast that structure formation 
via NEMPI is impossible. It is therefore unclear whether mean- 



6. Conclusions and discussion 

The present work has confirmed the rather stringent restric- 
tions of LBKMR showing that NEMPI is suppressed already 
for rather weak rotation (Co 0.03). Applied to the Sun with 
f2 = 2 x 10~ 6 s _1 , this implies that NEMPI can only play a role 
in the upper-most layers where the correlation time is shorter 
than Co/2f2 w 2 hours. This is surprising as it does exclude even 
the lower parts of the supergranulation layer where r is of the or- 
der of 1 day. On the other hand, we have to keep in mind that our 
conclusions based on isothermal models, should be taken with 
care. It would therefore be useful to extend the present studies 
to polytropic layers where the scale height varies with depth. It 
should also be noted that weak rotation (Co = 0.03) seems to 
enhance the surface appearance. This would be compatible with 
an interpretation that rotation prevents structures from sinking. 
However, preliminary mean-field calculations do not indicate a 
significant dependence of the eigenfunction on Co for values be- 
low 0.1. Our interpretation, if correct, would therefore need to 
be a result of nonlinearity. 

If we were to apply NEMPI to the formation of active re- 
gions in the Sun, we should keep in mind that the scale of struc- 
tures would be 6-8 pressure scale heights (Kernel et al., 2012a). 
At the depth where the turnover time is about 2 hours, we es- 
timate the rms velocity to be about 500 m/s, so the scale height 
would be about 3 Mm, corresponding to a NEMPI scale of at 
least 20 Mm. This might still be of interest for explaining plage 
regions in the Sun. Clearly, more work using realistic models 
would be required for making more conclusive statements. 

Regarding the production of kinetic helicity and the possible 
detection of a magnetic helicity spectrum, our results suggest 
that the relative magnetic helicity cannot be expected to be more 
than about 0.01. This is a consequence of correspondingly low 
values of kinetic helicity. We find that the normalized kinetic 
helicity is given by ef « 2GrCo. For the Sun, we expect Gr = 
(fcfiJp) -1 w 0.16 (see Appendix C), which agrees with what 
is used in our simulations, leaving therefore not much room for 
more optimistic estimates. 



Acknowledgements. We thank Koen Kernel for helpful comments concerning 
the influence of magnetic structure formation on the measurement of q p o and 
/3 P in DNS. Ilia R. Losada was supported by PhD Grant 'Beca de Investigation 
CajaCanarias para Postgraduados 2011'. This work was supported in part by the 
European Research Council under the AstroDyn Research Project No. 227952, 
by the National Science Foundation under Grant No. NSF PHY05-51 164 (AB), 
by EU COST Action MP0806, by the European Research Council under the 
Atmospheric Research Project No. 227915, and by a grant from the Government 
of the Russian Federation under contract No. 1 1.G34.3 1.0048 (NK, IR). We 
acknowledge the allocation of computing resources provided by the Swedish 
National Allocations Committee at the Center for Parallel Computers at the 
Royal Institute of Technology in Stockholm and the National Supercomputer 
Centers in Linkoping. 



8 



I. R. Losada et al.: Competing effects in concentrating magnetic flux 



0.001 



Co 



0.01 




C 2 = A 2 - 7A 1 + 35Ci , 
C 3 = A x - 5C\ . 

In the case of a <C 1 these functions are given by 



(B.4) 



Mo) 

Ci(o) 



47T 

T 

4-7T 

15 



Ma) 



8tt 
15 



a , 



0.01 



2fi/A 



*o 



In the case of a » 1 these functions are given by 

7T 2 87T 

y/a a 

CM 

The functions An (/3) are given by 



IT 2 


47T 


sja 


1 

a 


7T 2 


47T 


4^ 


3a 



3/3 6 fP Rml/i A n (X 2 ) 



X 7 



dX, 



2tt 



Fig. A.l. Dependence of A/A»o on 2S7/A*o for 9 = 0° in the 3D and 
2D cases for — (corresponding to the pole). 

Appendix A: Revised comparison of the 2-D and 
3-D data of LBKMR 

We note that in Fig. 3 of LBKMR the growth rates from the 2- 
D calculations came out to be much smaller than in 3-D. This 
was however a mistake caused by determining the amplification 
of the total magnetic field (which included the imposed field), 

rather than the fluctuations of the magnetic field or the veloc- plicit form of the functions A„ 0) and C„ (/?) for m — 0; 2 
ity. We use here the opportunity to present the revised plot in are given by 



! A n 
Jo 



1/4 



A n (X 2 ) 
X 3 



dX 



(a(f))f m df= Y A^0), 



(B.5) 
(B.6) 
(B.7) 



and similarly for ci" l) (/3), where a = [f3u kT(k)/2] 2 , (3 = 
y/8 B/B cq , and X 2 = f3 2 {k/k ) 2 ^ = /3 2 /f = a. The ex- 



Fig. A.l. 

Appendix B: The identities used in Sect. 4.2 for 
the integration in k-space 

To integrate over the angles in fe-space in Sect. 4.2 we used the 
following identities (Rogachevskii & Kleeorin, 2004, 2007): 

ka sin 9 



K, 



1 + a 



cos 2 



dd dip = AiSij + A 2 pij, 



(B.l) 



K, 



kijmn sin 
1 + a cos 2 6 



d6 dip = C\(Si 3 S 7nn + Si m Sj n 



-\-Si n S jm) C2^ijmn "T~ C^^Sij [3 7nn -{- Simftjn ~\~ $in$jr 



1 1 j run 



■^i?'mn(^) ' 



2n,2 dddl P 



(B.2) 



kijmn sin ( 



(1 + acos : 

/kijmn sin 6 
b + a cos z 



4 0) C9) 



4 0) (/3) 



A ( ?\P) 



cf(£) 



arctan /3,„ 
2 + 2 ? — (3 + 5/3 2 ) 



G 



-2/3 2 ln 
2 " 

- 2- 



/?3 

1 + /3 2 
l + /3 2 VRrri 
arctan /3 



/3 2 InRm 



(B.8) 



(9 + 5/3 2 ) + -/3 2 lnRm 
/3 3 V ; /3 2 



2 

63 



2/3 2 In 
1 + 



1 



3 Z 



1 + /3 2 VRm 
arctan j3 
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where = ftj = /3i/3j, and 

2tt 



(B.3) 



a 

J 2tt 
A 2 = 



^ " 2V 



arctan(^)_ i 
v a 

(Q + 3) ar C tan(V5) _ 3 
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. _ , 9 arctanf-i/a) 5a 
(« + !) 



05 

9 /3 2 
2/I 2 ~ 2^"T 

where M(/3) = 1 - 2/3 2 + 2^ 4 ln(l + /3~ 2 ). Here we have taken 
into account that Rm >• 1. For £> <c i? cq /4Rm 1 / 4 these func- 
tions are given by 

4 0) (/3) 



2 - i/3 2 InRm, 
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For B oq /4Rm 1 / 4 <C B <C -B oq /4 these functions are given by 



4 0) (/3) ~ 2 + p 2 



15 7 



4° } (/3) - j^ 2 



41n^---3/3 2 



Other functions in this case have the same asymptotics as in the 
case of B <C _B cq /4Rm 1/ ' 4 . For B > 73 oq /4 these functions are 
given by 



Af\p) 
cf\p) 



/3 /3 2 2 V ' f3 2 



3tt 3 

7/3 + W 



3tt 3 
7/3 ~ 2^2 
3tt 1 

28/3 2/3 2 ' 
The functions q 1 (/3) and g 2 (^) are given by 

qi (J3) = 1 [4 0) (0) - ^i 0) (4/3) - i4 0) (4/3) 
-(Qt ) 2 (4 2) (0) - 2Cf } (0) - 104 2) (4/3) 
+206^(4/3) + A [^(ie/3 2 ) - 2C\(16/3 2 )] ) 



(B.12) 



= 4^4° ) (4/3) + (^ro) 2 (4 2) (0)-6Cf ) (0) 



12 



(B.13) 



-104 2) (4/3) + 60Cr J (4/3) + — L4!(16/3 2 ) 
-6Ci(16^ 2 ) ' 
where /3 = B/B eq . 

Appendix C: Revised estimate for k f H p 

In Kernel et al. (2012a) the value of k{H p (= Gr _1 ) was es- 
timated based on stellar mixing length theory, using £ m [ x = 
ctmixHp for the mixing length with a m i x H p w 1.6 as an em- 
pirical parameter. For isentropic stratification, the pressure scale 
height Hp is related to H p via jH p w i/ p . With fcf = 27r/£ m i x 
we obtain fcfifp = 27r7/a m i x w 2tt, so Gr = (kfHp)^ 1 ss 
0.16. Owing to a mistake, we underestimated in Kernel et al. 
(2012a) the value of kf H p by a factor of 2.6. The correct value of 
the growth rate of NEMPI should therefore be larger and would 
now be clearly faster than the turbulent-diffusive rate. 
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